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Abstract 

We present an analysis of transport properties of a system consisting of two half- 
infinite interacting one-dimensional wires connected to a single fermionic site, the 
energy of which is subject to a periodic time modulation. Using the properties of the 
exactly solvable Toulouse point we derive an integral equation for the localised level 
Keldysh Green's function which governs the behaviour of the linear conductance. 
We investigate this equation numerically and analytically in various limits. The 
period-averaged conductance G displays a surprisingly rich behaviour depending 
on the parameters of the system. The most prominent feature is the emergence of 
an intermediate temperature regime at low frequencies, where G is proportional to 
the line width of the respective static conductance saturating at a non-universal 
frequency dependent value at lower temperatures. 
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1 Introduction 



The future advances of the microelectronics depend crucially on the devel- 
opments in the field of miniaturisation of electronic circuits. The ultimate 
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basis devices are single molecules or even atoms supplemented by three elec- 
trodes: source, drain and an additional gating electrode responsible for the 
current flow control. Such structures can, however, be so small that quan- 
tum effects have to be taken into account. The simplest possible realisation 
of such a device is a single fermionic level coupled to two metallic electrodes 
and subjected to an electrostatic interaction with a gate electrode. Such struc- 
tures are also called single-state quantum dots. In the ideal situation of very 
small lateral device dimensions the contacting electrodes can be considered to 
be one-dimensional (ID). For instance, thin metallic single-wall carbon nan- 
otubes (SWNTs) can play the role of such wirings. However, the low energy 
sector of ID metals cannot be described by the conventional Fermi liquid the- 
ory, they are instead Luttinger liquids (LLs, for a recent review see e.g. [1]). 
Generally, the transport problem of a quantum dot coupled to LLs is not inte- 
grable [2]. Nevertheless, in some situations there exist unitary transformations 
which allow one to map the Hamiltonian of the system onto a diagonalisable 
quadratic one. In such situations all transport properties including the noise 
power spectra are directly accessible [3,4]. 

In addition to the static conductance properties future applications in inte- 
grated circuits demand also a detailed knowledge of the device response to 
AC currents and especially to an AC modulated gating. Thus far numerous 
studies of quantum dots under such conditions have been made [5,6,7]. How- 
ever, all of them deal with quantum dot systems contacted by non-interacting 
electrodes. The purpose of this work is to gain initial insights into phenomena 
taking place in AC-modulated devices coupled to interacting LL electrodes. 
The outline of the paper is as follows. We begin by summarising the results for 
non-interacting systems in Section 2, where we also elaborate on the mean- 
ing of the low and high frequency approximations. We analyse the dynamic 
conductance for LL electrodes in Section 3. 



2 Formulation of the problem and non-interacting results 

The Hamiltonian under investigation, 



contains three parts. The first one, Hk, describes the dynamics of the resonant 
level and the electrons in the right/left (R/L) electrodes (we shall also call 
them leads), 



H = H K + H t + H c , 




H K = A(t)cftd+ E HoN, 



i=R,L 
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where the dynamic gating of the resonant level is described by explicitly time 
dependent level energy A(t). As long as the Kondo temperature is small, the 
spin effects are trivial for resonant tunnelling problem, therefore we work with 
a spinless model and briefly discuss the results for the spinful situations in 
Conclusions. The tunnelling between the resonant level and the electrodes 
with energy independent amplitudes 7r 5 l is taken care of by the second term, 

if t = 5>i[<*tyi(0)+h.c.]. 

i 

In (1), He stands for the electrostatic Coulomb interaction between the leads 
and the dot, 

H c = \ c Sd ^1(0)^(0) . 

i 

In the non-interacting case, when H [if)i] describes free fermions in ID and 
when Ac = the non-linear I(V) can be shown in the static case to be given 
by the Meir-Wingreen formula [5], 

I(V) = — /^[/LH-^H]-^^ImL>>), (2) 

IT J L R + I L 

where /l,_r(u;) = n R {uJ z feV/2) are the Fermi distribution functions in the leads 
and D r (lo) denotes the Fourier component of the retarded Green's function 
D R (t, t') of the localised level. T RjL = 7rp 7i? l are the dimensionless conduc- 
tances of the contacts (p denotes the density of states in the leads). For the 
linear conductance one obtains then 

G/G = ^ J duo cosh' 2 (cu/2T)T(uj) (3) 

where T(u) = —T L T R lmD R (uj)/(T R + T L ) is the transmission coefficient of 
the electrons through the double constriction and Go = e 2 /h denotes the con- 
ductance quantum. In the static case D R (t,t') depends only on the difference 
t — t' while in case of the dynamic gating it depends on t — t' as well as t + t'. 
In such situations one is usually interested in the current (or conductance) 
averaged over the period of the gate voltage oscillations. The non-interacting 
situation has been discussed in detail in Ref. [6]. Here we are interested in 
the linear conductance rather then in the full I{V). Therefore we recapitulate 
and somewhat extend the results of Ref. [6] in order to be able to compare to 
the interacting results we're presenting in the next Section. In this situation 
D R (t, t') can be solved for exactly leading to [7] 
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Fig. 1. Ze/ii 'panel: Temperature dependence of the period averaged linear conduc- 
tance for different gate oscillation frequencies and A\jY = 5, Ao = 0. The thick 
dotted line (not included into the legend) represents the result of the adiabatic 
approximation. Right panel: Thick lines: frequency dependence of the asymptotic 
linear conductance at T = for different values of Ai, Ao = 0. Thin lines represent 
the cases of nonzero Aq/T = 2, 5, 8 and Ai/T = 5 (from above). Notice the double 
logarithmic scale. 

t 



D R {t,f) = -iQ(t-t')exp 



-i J drA(r)-^(T L + T R )(t-t') 



(4) 



Assuming the gating is harmonic, A(t) = A + Aicos(fit), one can easily 
perform the time integrations and the averaging over the period 27r/Q with 
respect to t + 1' . Then, using the Bessel Jq function one obtains [6]: 



D R (t - t') = -iQ{t - f) e [-<Ao-(r i+ r«)/2]( t -f) Jo j^i sin[fi(t _ t >y 2 ^ . 

As a result we obtain for the imaginary part [from now on we measure all 
energies in units of (T R + T L )/2] 



lmD R (u) = - £ 



J^AJQ) 



(lu-A + nil) 2 + 1 



(5) 



Insertion of this relation into (3) yields then the following analytic formula for 
the linear conductance 



G/G = 



^(A^Re *' 



1 1 + i(ntt - A ) 

2 + 2rfT 



•(6) 



At least in the non-interacting situation the asymmetry effects are trivial and 
are only reflected by the corresponding prefactor, so that we restrict our anal- 
ysis to = T^. At T ^> T, the conductance ceases to depend on either Q or 
Aj (i = 1, 2), and is given by G/G ~ 7r/4T. In the opposite limit of low tem- 
peratures the main feature of G is generic: it saturates as T approaches zero. 
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However, the saturation values are not the same at different oscillation fre- 
quencies. For very small frequencies and A < A 1; G(0)/Gq w 1/Ai whereas 
for A >> Ai, the conductance decays as ~ Aq 2 , see Fig.l. From the physical 
point of view it can be understood as follows. In the static situation due to 
finite tunnelling amplitude the localised level and the leads hybridise, so that 
the spectral function of the dot state has a finite width ~ T. Then, as long as 
Ao < Ai, the level spends a fraction v — V j A\ of the oscillation period within 
'reach' of the electrodes. That's why the averaged conductance is proportional 
to v. On the other hand, for A > A 1 the overlap of two areas cannot ever be- 
come maximal. That leads to a fast decay of the averaged G as a function of A . 
In the opposite case of very high frequencies the average time, which an elec- 
tron needs to be transferred through the constriction, 1/r, is much larger than 
the oscillation period, so that the particles only see the level being 'frozen' at 
its average energy A , which means that G/G « (1 + A 2 ,) -1 in that limit. The 
crossover between large and small frequencies is highly non-monotonous, see 
Fig.l. There are local maxima at Q = Ao/(2n + 1) and minima at Q = Ao/2n 
{n is a natural number), which correspond to absorption of multiple energy 
quanta. 

Of course, all these results can be recovered using the corresponding asymp- 
totic expansions of the exact result (6), see Appendix. Moreover, the asymp- 
totic behaviour of the averaged conductance can be extracted from the static 
level Green's function without knowledge of the exact dynamic one (4). At 
very low frequencies O C T the system can be regarded to be in equilib- 
rium with respect to the oscillations at any time and the obvious approxi- 
mation is thus given by (2) with the period-averaged transmission coefficient 
T(uj, t) — [1 + {uo — A(t)) 2 ] -1 (we are still dealing with the symmetric situation 

r R = r L ) 



2-k/Q, 

T(cu) = |- J dtT(uj,t) = Re{[l + t(uj-A )] 2 + Aiy 1/2 . (7) 
o 

The results of this adiabatic approximation are depicted in Fig.l and describes 
very well the behaviour of G for frequencies up to T. We wish to stress here 
that Q — > limit is non-trivial and is not equivalent to simply taking Ai = 0. 
Since the usual experimental values of T range between 0.1 and 1 /xeV, the 
static approximation is expected to be valid to fairly high frequencies of the 
order of several MHz. The opposite limit of Q ^> T is even simpler. In that 
situation one can use the formula for the static case with A(t) — > A . This 
approach we refer to as anti- adiabatic approximation. 
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3 Interacting leads 

It has been shown in a number of recent works that the resonant tunnelling 
between correlated electrodes is completely different from the non-interacting 
case [8,9,10,3,11,4,12]. As in the non-interacting situation, we start with the 
Hamiltonian (1) with Ac* 7^ 0. However, contrary to the previous case, we 
model the leads by half-infinite LLs. It turns out that in the static case the 
problem can be solved exactly at the so-called Toulouse point, when the LL 
interaction parameter g — 1/2 and Ac = 2n [3]. Generally, the transport 
properties can be extracted in two different ways, either using the equations of 
motion method [3], or via Keldysh diagram technique [4]. The latter approach 
is, however, considerably easier to apply in the case of a dynamic resonant 
level. After a unitary transformation, refermionization and introduction of 
the Majorana components of the new fermionic degrees of freedom, £, rj, and 
those for the dot level, a, b, the resulting Hamiltonian can be written in the 
form H = H + H t , where (see [3,4] for details of the mapping; we use the 
same notation) 



and H t = i7&£(0), while the current operator is given by J = —276^(0). 
Calculating the average of this operator using the S matrix containing the H t 
coupling, one can derive the following expression for the full (time-dependent) 
non-linear I — V of the system, 



where T = 7 2 and D bb (uo,i) is the Fourier component of the retarded Green's 
function (GF) of the b Majoranas with respect to the difference of time argu- 
ments, 



As in the non-interacting case, the linear conductance turns out to be given 
by the same formula (3). The meaning of the transmission coefficient T(u), 
of course, changes. Now it is related to the probability for the 'new' fermions 
(see Ref.[3]) to be scattered off the resonant level. According to Ref. [4] in the 
static case one obtains 







(10) 



TH=c 2 /[(c 2 -A 2 ) 2 + c 2 ], 



(11) 
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Fig. 2. Temperature dependence of the averaged conductance in the interacting case 
for different gate frequencies Q and Ao = 0. The uppermost curve represents the 
result of the anti-adiabatic approximation while the lowermost curve is the result 
of the adiabatic approximation. The curves in between correspond to Q/T = 4, 3.5, 
3, 1, 0.1 and 0.05 (from above). Notice the double logarithmic scale. 

while for dynamic gating the effective T(oj) has to be found via solution of 
corresponding Dyson equation for the retarded GF, 



D&t, t') = £>£>*(*, t') + Y J drdr'D^ R (t, t)G$\t - t')D*(t, t') , (12) 



where (r — r') = — i6(r — t')/2 is the retarded GF of the lead Majoranas 
for 7 = and 



is the bare retarded GF of the resonant level Majoranas. Unfortunately we 
could not solve Eq.(12) analytically. Before we present numerical results, let 
us discuss the limiting cases of low and high Q. As expected from the analysis 
of the non-interacting case, these limiting cases should be reliably described 
by the adiabatic and the anti-adiabatic approximations, respectively. Let us 
first investigate the most interesting case of zero off-set A = 0. The f2 — > oo 
limiting case is trivial and the conductance is given by the corresponding static 
result with Ao,i = 0, see the uppermost curve in Fig.2. In the case of slow os- 
cillations we resort to the adiabatic approximation, which can be implemented 
in the same way as in the non-interacting situation, namely by averaging the 
corresponding static transmission coefficient (11) over the oscillation period: 




(13) 



t> 




2u 



Im{[2(cu 2 



)/A?-l] 2 -l} 



1/2 



(14) 




— 10J 



Evaluation of the linear conductance using this approximation results in the 
lowest curve in Fig. 2. The main feature is that G vanishes at zero temperature 



as ~ yT/r. To explain that we can proceed along the same line of reasoning 
as in the non-interacting case. It is reasonable to assume that in this particu- 
lar situation the dominant contribution to G is generated near the resonance, 
when A is aligned with the chemical potentials in the leads. The averaged 
conductance of the structure is then proportional to the ratio of the effective 
level width w(T) to the amplitude A 1 of the oscillation G ~ w(T)/Ai. As 
we know from the solution of the symmetric static case at g — 1/2, the level 
width vanishes as a square root of temperature [10,3]. That is why we find 
the same power-law for G(T). Note that this is very much different from the 
non-interacting case where, of course, w(T) does not depend on temperature. 
In order to obtain information about G at intermediate frequencies one has 
to solve equation (12) and extract the transmission coefficient numerically. 
The outcome of this approach is shown in Fig. 2. As can be seen from this 
figure, there is an important qualitative correction to the adiabatic approx- 
imation. Namely, the exact G saturates as T — > at a finite ^-dependent 
value: G(T = 0) ~ ^/u. The on-set of this saturation takes place at T* ~ Q. 
At temperatures higher than T* the system finds itself in the low frequency 
limit and G follows the square root power-law found in the static approxima- 
tion. On the contrary, below T* the constriction can be considered as being in 
the high frequency limit, where the conductance is expected to saturate. The 
resonant constellation A(t) ~ ceases to generate the dominant contribution 
becoming too narrow in comparison to the off-resonant positions of the dot 
level, when the effective width is known to saturate [10,3]. 



The results for the 'immersed' level, when Ao < Ai are very similar, see Fig. 3. 
Again, the conductance curves saturate at temperatures T <C Q. However, the 
(anti-)adiabatic approximation as well as the dynamic curves in the intermedi- 
ate regime show different power-laws: instead the square root in the resonant 
(Ao = 0) situation the conductance decays toward low temperatures quadrati- 
cally. This power-law originates in the different transport mechanism: now the 
most of the conductance does not come from the resonant transmission, which 
cannot be achieved any more, but from the off-resonance processes which show 
quadratic temperature behaviour. Moreover, in the immersed case the adia- 
batic approximation yields higher conductances than the anti-adiabatic ones. 
The asymptotic conductance behaviour at T = shows another interesting 
feature. It is maximal when the oscillation frequency matches the off-set energy 
A . In that situation the tunnelling electron can absorb an energy quantum 
from the gate and get transmitted through the system as if there is no energy 
off-set. 
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Fig. 3. Temperature dependence of the averaged conductance in the interacting case 
for different gate frequencies f2 and Ao/r = 2., Ai/T = 1. . The upper solid line 
represents the result of the adiabatic approximation while the lower one corresponds 
to the anti-adiabatic approximation. Notice the double logarithmic scale. 

4 Conclusions 



We have analysed transport properties of a single site quantum dot coupled to 
two LL electrodes in the case of an AC modulated gate. Using a unitary trans- 
formation and refermionization we were able to recast the Hamiltonian into a 
quadratic form at a special interaction strength g — 1/2. This allowed for an 
exact calculation of the experimentally relevant period averaged conductance. 
Its temperature dependence turns out to possess a rich behaviour, showing 
two cross-over points. Upon lowering the temperature for T 3> V the effect of 
the AC modulation is not visible and G(T) follows the same power-law as in 
the static case ~ T/T no matter how large is the off-set energy A . In the case 
of A = and for oscillation frequencies Q ^> T the averaged conductance 
saturates at a value larger than G(T = T). However, as soon as Q becomes 
smaller than T, the conductance vanishes first as a square root of temperature 
in the interval between T = f2 and T = T and after that saturates at a finite 
value. The saturation toward low temperature is the most prominent result 
of a dynamic gating and has a potential to obscure the resonant tunnelling 
physics occurring in strictly static set-ups. Nevertheless, an AC modulation 
experiment has additional advantages against the static one - it enables an ex- 
traction of the temperature dependence of the resonance peak width from the 
information about G(T) in the intermediate temperature regime VL < T < T. 

Our findings are, strictly speaking, restricted to the spinless g = 1/2 case. 
Based on this results and the achieved understanding of various conductance 
mechanisms, we can however speculate about the power-laws that will emerge 
in the general-^ case as well as in set-ups involving electrons with spin (so 
long as Kondo physics is not important) and SWNTs. For the spinless system 
the most important regimes are the following: (i) the high temperature limit 
G{T) ~ T 2g ~ 2 ; (ii) the intermediate regime G(T) ~ T 1 ^ 9 ; (iii) saturation 
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to a frequency dependent value at lowest temperatures. The results for the 
resonant tunnelling between two spinful LLs and SWNTs are obtained by a 
substitution g — > (g + l)/2 for a system with spin and g — > (K + 3)/4 for an 
SWNT, where K is the nanotube Luttinger parameter. We shall present the 
detailed analysis of the general-^ case in a separate publication [13]. 

To summarise, the dynamic resonant conductance for LL leads is markedly dif- 
ferent from that for non-interacting contacts. It possesses various distinctive 
features which are, in principle, verifiable experimentally. In our view a com- 
parison of static and dynamic (AC-modulated) conductance measurements on 
the same system would be especially valuable. 



A Adiabatic approximation 



In this Appendix, we elaborate on the exact formula (6) for the conductance 
in the non-interacting case at T = 0. Expanding (6) in T we obtain: 



First of all, it is useful to notice that G/Go, as given by formula (A.l), is 
always smaller than 1. This immediately follows from the well known identity 
(see, e.g., [14]) 

oo 

E ■£(*) = !■ 



The anti-adiabatic case (Q ^> A x ) is trivial. For A = the perfect con- 
ductance is recovered. Expanding the Bessel functions, we find the next-to- 
leading correction as 



Note that the perturbative correction in Ai is converging for all fl. 

Next we consider the adiabatic case z = fi/Ai <C 1. This turns out to be more 
involved. Using another well known identity (see [14]), 

oo 

E fi(z)e- in * = Jo[2^sin(0/2)], 

n=— oo 
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we re-write formula (A.l) as 

G/G = -Re / #J o [2zsin(0/2)]/(0) , (A.2) 

7T J 



with 

D ind> 



e" 



(ntt - A ) 2 + 1 



This sum can be evaluated by using Poisson summation formula with the 
result 

n e (7r-0)(l-iAo)/n n e (^.-7r)(l+iAo)/n 

^ = fi2sinh[7r(l -iA Q )/fl] + Q 2 sinh[7r(l + iA )/Q] ' 

When substituting the above expression into formula (A.2) for conductance 
we observe that, in the Q — > limit, the ^-integral converges in a narrow 
interval ~ Q. Rescaling <fi = Q9 and taking Q — > limit we find 

oo 

lim G/G = Re / d0 e - (1 - iA ° )fl J (A^) = Re . . (A.3) 

I ^(l-iA ) 2 + A? 



This result justifies the use of the adiabatic approximation adopted in the 
main text. 
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